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ABSTRACT 

We present a measurement of the velocity flow of the local universe relative to the CMB rest frame, based on the 
Jha, Riess & Kirshner (2007) sample of 133 low redshift type la supernovae. At a depth of 4500 km s -1 we find 
a dipole amplitude of 279 ± 68km s" 1 in the direction / = 285° ± 18°, b = -10° ± 15°, consistent with earlier 
measurements and with the assumption that the local velocity field is dominated by the Great Attractor region. 
At a larger depth of 5900 km s -1 we find a shift in the dipole direction towards the Shapley concentration. We 
also present the first measurement of the quadrupole term in the local velocity flow at these depths. Finally, we 
have performed detailed studies based on N-body simulations of the expected precision with which the lowest 
multipoles in the velocity field can be measured out to redshifts of order 0.1. Our mock catalogues are in good 
agreement with current observations, and demonstrate that our results are robust with respect to assumptions 
about the influence of local environment on the type la supernova rate. 
Subject headings: 



1. INTRODUCTION 

Distant type la supernovae (SNe la) have been used to probe 
the expan sion rate of the U niverse out to redshifts of order 1.5 
(see e.g. lRiess et al.ll2004b . These measurements were crucial 
in establishing the current standard model of cosmology in 
which roughly 30% of the energy density is in the form of 
non-relativistic matter, whereas roughly 70% is in the form of 
a dark energy, a component with negative pressure. 

The power of supernova surveys as cosmological probes 
depends on a precise measurement of the luminosity dis- 
tance to the individual supernovae - but the supernova host 
galaxies do not follow the Hubble flow. They have peculiar 
velocities, induced by the underlying gravitational potential, 
and the measured luminosity distances and redshifts are per- 
turbed in correlation w ith the large scale structure (see e.g. 
iMiller & Brar^fl992h . 

The measured fluctuations in the luminosity distances from 
nearby supernovae surveys ca n be related to the local varia- 
tion in the Hubble parameter (Riess et al . 1995), and the lo- 
cal large scale structure, and can be quantified in terms of 
their correlation functio ns. This possibility h as been studied 
in se veral recen t paper s ( Bonvin et al. 2006a ; b; Hui & Greene 
2006) (see also[Su siura, Susivama & Sasaki 1999) using an- 
alytical methods. iBonvin et al.l (|2006b|) detected the dipole 
term with respect to the CMB at roughly 2a confidence. The 
variati on of the monopole c ontribution with redshift was stud- 
ied by Zehavi et~aD (1 1998b using the same supernova sample 
as (|Riess et al.lll995|) . 

In contrast to measurements of the density field, velocity 
field measurements are much less sensitive to selection bias, 
since the field is measured directly instead of summed up from 
number counting, but it is very sensitive to uncertainties in the 
measured luminosity distance. Because of this type la super- 
novae are particularly useful as probes of the velocity field 
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because of the very small inherent uncertainty in their lumi- 
nosity distance. Compared with measurements using galaxies 
as standard candles much fewer supernovae are needed to get 
a reliable estimate of at least the lowest multipoles in the ve- 
locity field. 

The purpose of the present paper is twofold. The first is to 
perform a theoretical study of the problem of extracting the 
velocity field from supernova data. The second is to use the 
developed formalism on the best available data set to extract 
precise values of the dipole and quadrupole terms in the local 
velocity field. 

We first use large scale dark matter N-body models to pre- 
dict the observed angular power spectrum of the peculiar ra- 
dial velocity field as a function of redshift and to explore its 
utility to probe the local velocity field. By using this method 
we are able to quantify the effect of various error sources as 
well as cosmic varia nce. 

We also use the iJha. Riess~& Kirshner (2007) (hereafter 
JRK) sample of nearby type la supernovae to calculate the 
lowest multipoles of the local velocity field. We find that both 
the dipole and the quadrupole terms are well measured by 
this sample. Higher multipoles cannot be reliably estimated 
with t he JRK sample because of sparse sampling. In I Jha et all 
( 2007) this sample was also use d to probe the mono pole term 
in an analysis similar to that of IZehavi et al.l (fl998|). Similar 
evidence of a local void was found in Jha, Riess & Kirshner 
(|2007|) . even taking a data set w hich is completely disjoint 
from that in IZehavi et al.l {1998). In the following we will 
concentrate on the dipole and quadrupole terms in the JRK 
sample si nce the monopol e has already been exhaustively dis- 
cussed in Jha et al.l (200% 

In Section 2 we discuss the formalism used to derive the 
velocity field from magnitude measurements. In Section 3 we 
present the analysis tools used to extract the multipole compo- 
nents of the velocity field from a finite sample of supernovae. 
In Section 4 we use catalogues from N-body simulations to 
study the expected properties of the velocity field, including 
the precision with which it can be probed by supernova sur- 
veys. We apply the same formalism to the JRK sample and 
provide a detailed discussion of the results in Section 5. Fi- 
nally, Section 6 contains a comparison with other measure- 
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ments of the local velocity field, and in Section 7 we provide 
our conclusions. 
In the following we will take c = 1 . 

2. DISTANCE MEASURES AND RADIAL VELOCITIES 

The luminosity distance, d^ to a supernova at redshift z is 
given by the definition 



d r = 



4ttF ' 



(1) 



where F is the observed flux and L the luminosity, or equiva- 
lently by 

w=5i ° g -(^y +M+25 ' (2) 

where m is the apparent and M the absolute magnitude. The 
angular-diameter distance, d A , to the same supernova is de- 
fined as 

In a homogeneous and isotropic universe the two distance 
measures, and d A , are given by 



d A (l+z) = 



1 



(l+z) H 



where 



H [ 

Jo 



z dz 
o H(z) 



l+z 



Q>1 

n = i , 

sinh(^Vl-^/) fi<l 

(4) 



(5) 



By taking the logarithmic derivative of Eq. © with respect 
to (1 +z) in the special case of a flat universe = 1) we obtain 

adz)- 1 = a A (z) + 1 = ttt^tvt = xi^fl ; > ( 6 ) 



H(z)d A (z) H(z)d L (z) ' 



where 



dlndL 



and 



dlnd* 



(7) 



Jln(l+z) dln(l+z) 

However, in a perturbed universe, the luminosity distance 
depends on the detailed trajectory of the individual photons 
from the supernova to the observer. At moderately high red- 
shifts (z > 0.5) the contribution arising from lensing by in- 
tervening matter dominates, while at low redshifts (z < 0.5) 
the contribution from the peculiar velocities of the supernova 
host galaxy relative to the observer dominates. In this paper 
we are only concerned with the local universe. From now 
on we will therefore only consider the contribution from the 
peculiar velocities in determining distances. For a given su- 
pernova and observer, each with some peculiar velocity, the 
measured redshift z, angular-diameter distance, d^ , and lumi- 
nosity distance, dj^ are mo dified according to dBonvin et al.l 
l200^lHuT & Greendl2006h 

l+z=(l-v -fl)(l+z)(l+v r ) (8) 

d A = d A (z)(l+v -n) (9) 

d L = d L (z)(\+vo • n)(l + v r - v • nf , (10) 
where by definition 

d L = d A {\+zf (11) 

d L (z) = d A (z){\+zf (12) 



A bar .7. indicates the quantities as measured in a homoge- 
neous and isotropic cosmology, and v r = v e • n is the velocity 
of the supernova projected along the direction from the ob- 
server to the supernova. Having measured the redshift, z, and 
the flux, F, of the supernova we can calculate the luminosity 
distance, d^ and the angular-diameter distance, d A , to within 
a scatter determined by the cosmic variance on the luminos- 
ity, L, of a supernova. As we know the peculiar velocity, vo, 
of the observer with respect to the CMB with great accuracy 
it is useful to collect all known quantities on the left hand side 
of the three equations: 

l+z=(l+z)(l+v -/i) = (l+D(l+v r ) (13) 

d A = d A (l -v -n) = d A (z) (14) 

d L = d L (l + v • n) = d L (z)(\+2v r ) (15) 

The three quantities z, d A , and are the redshift, angular- 
diameter distance, and luminosity distance as measured and 
calculated by an observer at rest with respect to the CMB. In 
the following we assume that the measured quantities all have 
been corrected to a frame at rest with respect to the CMB. The 
measured redshift, z, is most transparently split into the cos- 
mological redshift, z, and the radial velocity of the supernova, 
v r , by inverting Eq. (ITU) to find z, which is then inserted into 
this alternative form of Eq. (fT3l) : 

v r = ln(l+z)-ln(l+z) (16) 

It is also possible to begin with Eq. (IT6b to obtain the cos- 
mological redshift, z, in terms of the measured redshift, z, and 
the radial velocity, v r , of the supernova: 

ln(l+z) = m(l+z)-v r (17) 

We now expand the natural logarithm of the angular-diameter 
distance, In d A (z), to first order around the measured redshift, 
z, and by using Eq. (ITU) we obtain 

lnJ A = lnd A (z) = lnd A (z)-a A (z)v r , (18) 

where a A (z) for a flat universe is given by Eq. ©. This equa- 
tion applies equally well for the luminosity distance, so we 
can give a common formula for both distance measures: 

In f J^r U In [ 4^r I = ~a A (z)v r (19) 



, d L (z) 



d A (z) 



This expression can easily be transformed into an equation 
between the measured apparent magnitude, m, the calculated 
apparent magnitude, m(z), at the measured redshift, z, and the 
peculiar radial velocity, v r , of the supernova: 



m - m(z) = 5 log 



10 



d L (z) 



= -(5/lnl0)a A (£)v r (20) 



Inverting this equation for the special case of a flat universe 
we obtain 

v r = -(lnl0/5) (J^§L\ (m-m(z)) (21) 
\l-H(z)d A (z)J 

For low redshift supernovae the Hubble parameter, the 
angular-diameter distance, and the luminosity distance can all 
be expanded in terms of the deceleration parameter, go : 

H(z)=H (l + (l+qo)z) (22) 



d A (z) = — (l-(3 + 90 )z/2) 
no 

d L (z)=^-(l + (l-qo)z/2) 
no 



(23) 
(24) 
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When we insert these values into Eq. (l2Tb we obtain the radial 
velocity of a low redshift supernova 

v r = -(In 10/5)z (1 + (1 + qo)z/2) (m-m(z)) (25) 

In this approximation the distance modulus, m-M, is given 
by 

m(z) -M= 42.3841 - 5 log 10 h 



+51og 10 z+(2.5/lnlO)(l-<7o)z, 



(26) 

where the Hubble constant as usual is given as Ho = 
100/zkm s _1 Mpc _1 . The distance modulus defined in Eq. (|26l) 
should be compared with the distance modulus, m-M, as 
measured in a frame at rest with respect to the CMB. The dif- 
ference between the apparent magnitude m and the one mea- 
sured by the observer, m, is according to Eq. ([T5K 

m-m- (5/mlO)vo -ft = 2.17vo -ft. (27) 

The approximate equation (l25t is sufficient for the low red- 
shifts that we are considering, while at higher redshifts one 
has to use the correct form Eq. (l2Tt . and also take into ac- 
count other contributions to m-m(z), such as lensing. 

3. ANALYSIS USING AN ANGULAR EXPANSION OF THE RADIAL 
VELOCITY FIELD 

We analyse both the mock catalogues discussed in the next 
section and the real data set (see section 0) using the same 
technique. In practice we decompose the field into spherical 
harmonics. 

The radial velocity is a real scalar field, and on a spherical 
shell of a given redshift it can be decomposed into spherical 
harmonics 



V r = ^ ^ aimYln 



/=() ra=-/ 



- ^ < -mXl -m + CllmYlm) + fl/pF/p 



(28) 



/=o 



Using ai - m = (-l) m aj m for the expansion of a real function and 
Y lrm = (-l) m Y? m we obtain 



CO ( I 



Vr = Y \ Y i^imYim)] + a m Y m \ (29) 

/=p I m = l J 

= £ ( Y l^imMYim) ~ 2Z(a lm mY lm )] + a/pF/p 1 

/=P t m = l J 

However, this applies strictly only if the field can be 
measured on the entire sphere. In our case the radial 
velocity field is measured for a finite number of direc- 
tions, so we can only hope to determine a finite number 
of ai m coefficients by fitting a truncated multipole expan- 
sion by the method of weighted linear least squares us- 
ing |Y/ , {23?(7/ m ),-2^(F /m )},m = 1, ...,/] as basis functions. 
Specifically, we solved the problem by a singular value de- 
composition. 

We follow the proced ure by 

ICopi. Huterer. Schwarz. and Starkmanl (|2006|) and repre- 
sent the /'th multipole in terms of a scalar, A (/) and / unit 
vectors, {v (/,m) ,m = 1, . . . ,/} : 



f l (0^)=A^{Y[(v^-e)-T l }, 



(30) 



where e = (sin#cos</>,sin#sin</>,cos60, and 7/ is the sum of 
all possible traces of the first term. In this representation the 
multipole expansion up to and including the quadrupole term 
takes the following form: 

v r ( e ~)=A (0) +A<V U) -e) 

+A( 2 >{(v (2 ' 1) - e ')(v (2 ' 2) -e)-(l/3)(v (2 ' 1) -v (2 ' 2) )} (31) 

Note that v (21) and v (2 2) are "headless" vectors only defining 
a line, not a direction. Equivalently they define a plane, but 
they do not define a rotation in that plane, so the normal to 
the plane is also headless. By convention we choose as the 
first vector, e\, the one with the largest absolute z-coordinate. 
We can choose e\ to point to the hemisphere near the pole 
without introducing a negative amplitude A (2) , if both e\ and 
£2 have their sign changed. Finally we define the normal to 
the plane spanned by the two vectors as ^ x This is the 
polar quadrupole vector. 

From the a\ m coefficients the monopole amplitude can be 
found as 

A «» = (32) 

V47T 



and the dipole amplitude and direction can be found as 

A (l) = (a 2 w + 2\a u \ 2 ) l/2 (33) 



> = -tan 



3?(an) 



-1 /«io \ 

:COS ( — J 



(34) 



(35) 



This is the direction of the maximum of the dipole. All the 
higher order multi pole vectors are found by using the program 
mpd_decomp by lCopi et alJ d2006h . 

For the lowest /-values the amplitudes in the multi pole vec- 
tor ex pansion are related to the usual power Q as (ICopi et alJ 
2006) 



Cb=47r[A (0) ] 2 



,1(2,1) . tf2,2)\ 



(36) 
(37) 

(38) 
(39) 

It should be noted that in general the individual multipole 
coefficients obtained in the fit to data can be strongly depen- 
dent on the number of modes included. The reason for this is 
that the window function does not cover the entire sky, rather 
there are patches with zero coverage. This means that the 
spherical harmonics are no longer orthogonal, and shows up 
as a leakage of power between different /. In fact, this is pre- 
dicted to be a significant problem for any harmonic analysis 
with limited sampling because the higher order multipoles do 
contribute significantly to the rms velocity. In Section [5TT1 we 
discuss the implications of sampling for the JRK sample. 

4. SYNTHETIC SUPERNOVA SURVEYS FROM MONTE CARLO 
SIMULATIONS 

Before analysing existing data we make mock catalogues 
of supernova data based on dark matter N-body simula- 
tions. This is done in order to get an estimate of the var- 
ious sources of error in such measurements. Th e N-body 
simul ations were done using the Gadg et- 2 code dSpringell 
120051: ISpringel. Yoshida & White! 12000) with a box size of 
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Peculiar velocities at e = 0.01 




Peculiar velocities at e = 0.03 




Peculiar velocities at e = 0.04 



-1726 _ " 1120 km 0"' 

Fig. 1 . — The variation in the peculiar velocity at z = 0.01 -0.04. 



800 Mpc h _1 and 512 3 & 768 3 particles respectively to make 
synthetic realisations. The box size is chosen large enough 
that the periodic nature of the box does not impact the simula- 
tion at the scales (z < 0. 1 , or < 300 Mpc h _1 ) we are interested 
in, and the high resolution run is made to assure that our re- 
sults are not dependent on the numerical resolution. In Fig.Q] 
is shown a typical all sky map in Mollweide projection of the 
peculiar velocity field at redshifts z = 0.01-0.04, computed in 
the CMB rest frame. 
The formation rate of SNe la as a function of the envi- 



ronment is not well known, though there are indications that 
at low redshifts the rate is directly propor tional to the stel- 
lar mass, a nd insensitive to the metalicitv (jNeill et al .1 12006) 
(though see ISharon et al.ll2006l: ISullivan et'al.ll2006l and ref- 
erences therein for indications of a bimodal distribution). On 
the ot her hand the semi- analytic estimates in th e literature (see 
e.g. iBonvin et al.ll2QQ6bt [gui & Greene 2006) assume a rate 
which is uniform on the sky. This is clearly not realistic, but 
to test the effect on the luminosity distance distributions we 
have made synthetic data sets using both a rate proportional 
to the mass and an uniform distribution. In Fig. [2] we show 
the distribution of peculiar radial velocities in the two mod- 
els. The differences between the two scenarios are at a percent 
level. Looking at Fig. Q] we see that the radial velocity field 
is smooth across voids in contrast to e.g. scalar fields like the 
density field. This is because it is only a pseudo-scalar, and 
the underlying vector field can be transported efficiently (i.e. it 
is easier to change the direction of a vector, than transport a 
scalar quantity). 

4.1. Making a mock supernova survey 

Real measurements rely on a tracer of the matter distribu- 
tion, whether it is galaxies or SN la. In any case, only a fi- 
nite number of objects will be available. Furthermore there 
will be a selection bias coming from the presence of galactic 
foregrounds etc. In this work we generate a mock supernova 
survey using the following strategy: 

a) The total number of measured SNe in a redshift bin, N, 
is chosen. 

b) Each SN is generated by sampling one of the two proba- 
bility distributions described in the previous section. 

c) To each SN a noise component is added, resulting from 
scatter in the (stretch corrected) intrinsic luminosity, uncer- 
tainty in extinction correction, measurement errors etc. We 
describe the errors as a Gaussian error with a spread of Am 
on the measured apparent magnitude of the SN. 

From this data set the angular power spectrum is calculated. 
For each type of simulated survey this task is performed 500 
times for 27 different observers, to find the mean and variance 
of the angular power spectrum. 

We choose a set ofN= 100 SNe per bin in 16 redshift bins 
at redshifts of 0.005-0.08 or equivalently with Hubble flow 
velocities 1500- 24000 km s -1 . This conforms roughly to 
the e xpectations from local supernova searches con ducted to- 
dav dJha et al.l200d:lKrisciunas et al.l2004l:lLi et al.l2003l) and 
in the near future ( Aldering et al.ll2002l: iFrieman et al.l l 2004; 
Hamuy et al. 2006), if they are binned into 3 or 4 redshift bins. 
In Figs.[3]|4]we show the evolution of the lowest multipoles as 
a function of z for both models of the supernova distribution. 
The red line and error bars show how a hypothetical survey 
without any external error sources, Am = 0, would perform. 
Hence, the error here is only due to the finite number of SNe 
that are used to probe the velocity field. The green line shows 
the same, but including a Gaussian scatter of Am = 0.08 is 
included. Using the 27 different realisations we can estimate 
the size of cosmic variance (shaded blue area). 

From the figures it is clear that with 100 homogeneous 
distributed SNe per redshift bin both the dipole and the 
quadrupole can be measured out to a redshift of about 0.1. 
Furthermore for the synthetic observations, we know the un- 
derlying cosmology and the real Hubble constant, and we can 
determine the monopole. In the case of real observations we 
can only measure the relative change. In other words the zero 
point is in principle only measurable asymptotically at high 
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Fig. 2. — The peculiar velocity distribution of supernovae for different 
types of scenarios, and at different redshifts. Sampling with a SNe la rate 
proportional to the density compared to uniformly on the sky only bias the 
velocity distribution a few percent. 



redshifts. Looking at the red curve, we see that for all mul- 
tipoles the power goes to zero at higher redshift, as the flow 
approach the background Hubble flow. We also note that cos- 
mic variance is very large in the monopole term. This in- 
dicates that th e local "Hubble Bubble" phenomena , such as 
those found in lJha et al.l (l2QQ7h : IZehavi et all d!998b . are not 
unlikely. 

The velocity amplitudes are positive definite, and includ- 
ing a Gaussian scatter in the velocities, is adding uncorrelated 
noise to the (synthetic) observations. Therefore the ampli- 
tudes of the multipoles, when errors are included, are overes- 
timated. This error is per se hard to separate from the signal, 
but can trivially be beaten down by a better control of the in- 
trinsic errors or increasing the number of SNe per redshift bin. 



Fig. 3. — Uniform distribution of Supernovae: The amplitudes of the mul- 
tipole vectors and the different errors in a SN la survey with 100 SNe per 
redshift bin and an intrinsic scatter in the magnitude of Am = 0.08. 



Given an observational data set synthetic observations with a 
realistic sky distribution should be used to separate the noise 
amplitude from the underlying velocity field. 

In accordance with the underlying velocity distribution 
(Fig. [2]), on average there is a 5% overestimation of the radial 
velocity amplitudes, when assuming the Sn la rate is propor- 
tional to mass compared to a uniform distribution. 

5. RESULTS FROM THE JRK SAMPLE OF NEARBY SUPERNOVAE 

We apply the analysis technique described in section [3] to 
the sample of 133 nearby supernovae obtained by JRK. The 
JRK sample includes 95 type la supernovae in the Hubble 
flow, with an intrinsic dispersion of less than 7% in distances. 
This is the best sample available with distances derived in 
a homogeneous way, by using the Multi color Ligh t Curv e 
Shape (MLCS2k2) method described in Una et al.l (120071) . 
This method is one of several others used to standardise su- 
pernova distances, but is particularly powerful in that it allows 
a disentangled correction for host galaxy extinction and also 
provides a statistically reliable way to estimate the errors. 

The entire JRK sample comprises 133 supernovae. We have 
selected 3 sub samples of these for our analysis. We followed 
iJha et alJ (120071) in selecting as the first sub sample 95 Hubble 
flow supernovae, the HF sample, a selection based on distance 
cut, requirement of an acceptable light curve fit and exclud- 
ing objects with very high extinction. The SNe in this sample 
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Fig. 4. — Supernova rate proportional to density: The amplitudes of the 
multipole vectors and the different errors in a SN la survey with 100 SNe per 
redshift bin and an intrinsic scatter in the magnitude of Am = 0.08. 



have redshifts between 0.0085 and 0.021 and a weighted av- 
erage redshift of z = 0.0196 or 5900 km s" 1 . The second sub 
sample, the 4500 sample, includes 74 SNe, and is similar to 
the HF sample, but without the highest redshift SNe and in- 
cluding a few SNe with a lower redshift. It has an weighted 
average of z = 0.015 or 4521 km s -1 and contains SNe with 
redshift between 0.007 and 0.035. The last sub sample, the 
3500 sample, includes 42 SNe, and is similar to the 4500 sam- 
ple, but without the highest redshift SNe. It has an weighted 
average of z = 0.0118 or 3550 km s -1 and contains SNe with 
redshift between 0.007 and 0.017. 

To the uncertainties in the distance moduli, we also add 
in quadrature an additional error of 0.08 mag in order to 
properly represent the final uncertainties following iJha et al.l 
d2007h . In addition we add 50 km s 1 in quadrature to the er- 
rors in the radial velocity in order to take the velocity disper- 
sion around the local ani sotropic Hubble flow into account. 
iKarachentsev et al.l (120031) find that the radial velocity disper- 
sion around the local (anisotropic) Hubble flow within 5 Mpc 
amounts to only 41kms _1 , when distance errors are taken 
into account. F or further details on the samples we refer to 
Uha etal.ld2007h . 
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Fig. 5. — The local universe as seen by the 2MASS sur- 
vey. The galaxy distribution image is courtesy of Dr. T.H. Jarrett 
(IPAC/Caltech) and the 2MASS team. Complete image can be found at: 
http://spider.ipac.caltech.edu/staff/jarrett/papers/LSS The arrows indicate 
important super clusters. 



In Fig. [9] we show the distance to the nearest supernova on 
the sky, measured in degrees, for all points on the sphere. 

The sample has good coverage, except for a few "holes" 
mainly defined by the galactic disk. The mean distance to the 
nearest supernova is 12.8°, but the maximum distance is 40.5° 
at b = 12.7°, / = 230.4°. There are three areas with distance 
larger than 30°, roughly centred on (b,l) = (0 o ,0°),(5°,80 o ), 
and (10°, 230°) respectively. For a given /, the distance be- 
tween zero points in the field is 180°//. If the largest holes 
in the sample have a size of A6 then the multipole decom- 
position becomes problematic around / ~ 180°/A#. For 
the present sample this corresponds to / ~ 2.5 so that the 
quadrupole can be robustly fitted, but not the octupole. We 
have tested this in practice. When the quadrupole is added, 
the monopole and dipole amplitudes and the dipole direction 
hardly change. However, when the octupole is included there 
is serious leakage of power from the lower /'s to / = 3 and 
the results change substantially. The reason is that the fit in 
the well-sampled regions can be improved by adding the ad- 
ditional 7 a^ m coefficients to the fit, but that this happens at 
the expense of very large changes in the unsampled regions. 
In order to probe the higher order multipoles it is essential to 
reduce the size of the voids in the sample. For a uniformly 
distributed sample of 95 supernovae the average distance to 
the nearest supernova is 10.4 degrees, nearly the same as in 
the JRK sample. However, the average maximum distance is 
30.1 degrees, significantly lower than in our sample. A uni- 
form sample of 95 supernovae could be used to probe / = 3 
robustly. For a uniform sample the average distance to the 
nearest supernova scales as AT 1 / 2 , and the average maximum 
distance in the sample roughly as AT - 4 . 

5.2. Results 

In Figs. [6l|8] we show the obtained 68% and 95% contours 
for the direction of the dipole and quadrupole vectors in galac- 
tic coordinates. The corresponding best fit values with their 
formal 68% errors are shown in Tables [T] and [2 

This result can be compared to other velocity sur- 
veys based on galaxy samples. The recent review by 
iFeldman. Sarkar & Watkins I (l2006|) summarises these sur- 
veys. At an effective depth of 4000 km s -1 they find that 
the dipole amplitude is 330 ± 101 km s -1 in direction / = 
234° ±11° and/?= 12° ±9°. 



5.1. The effective window function 



5.2.1. Dipole 
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FIG. 6. — The dipole vector (top) and polar quadrupole vector (bot- 
tom) calculated from supernova data, with an weighted average velocity of 
3500 km s _1 . The ellipses show the one and two sigma errors 




Fig. 7. — The dipole vector (top) and polar quadrupole vector (bottom) 
calculated from Sn la data, with an weighted average velocity of 4500 km s _1 . 
The ellipses show the one and two sigma errors 



The 3500 and 4500 samples both show a dipole in a direc- 
tion compatible with the Great Attractor region at / ~ 300°, 
b ~ 0°. For the HF sample the direction shifts to slightly 
higher b, compatible with a shift in the motion towards the 
Shapley concentration (which lies at an average distance of 
14000 km s' MBardelli etal. 1119941) at slightly higher galactic 
latitude. Furthermore the amplitude of the dipole decreases, 
as is expected from the Monte Carlo simulations. 

5.2.2. Quadrupole 

For all the three sub samples we find a relatively large con- 
tribution from the quadrupole, showing that the local flow has 
a significant shear component. The result is consistent in mag- 




Fig. 8. — The dipole vector (top) and polar quadrupole vector (bottom) 
calculated from the JRK "Hubble flow" sample. The ellipses show the one 
and two sigma errors 




Fig. 9. — The shortest distance on the sky to a supernova in degrees using 
the "Hubble Flow" sample 



nitude with the expectation from the Monte Carlo simulations. 
From the figures it can also be seen that there is a change in 
quadrupole direction with redshift, and that the distributions 
for the 3500 and 4500 samples are bimodal. This could be 
because the quadrupole is pointing in different directions at 
the lower and higher end of the included redshift range of the 
SNe, and hence for the 3500 sample the bimodality is more 
apparent. Even though we quote formal 68% errors in Table 
2 the distribution is highly non-Gaussian and the error bars 
should be taken as indicative only. 

6. COMPARISON WITH OTHER RESULTS 
6.1. SNIa 

ThelTonrv etal.l(l2003h data set of 98 SNe was analysed by 
iHudson et al.l (|2004|) in order to find the local dipole. For the 
part of the sample with v r < 6000 km s -1 they found a dipole 
of v r = 376 ± 81 km s" 1 towards / = 285°, b = -14°. They 
do not quote error bars on this result, but for the part of the 
sample with v r > 6000km s _1 the stated errors are ±17° for / 
and ±13° for b. This result is completely compatible with our 
result for the dipole. However, Hudson et al. quote no results 
for the higher order terms. 
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Fig. 10. — The one dimensional distributions of the velocity amplitudes 
for the different samples, with the 68% limits indicated by a dashed line and 
the median value by a dotted line. 



Longitude for the 3500 km s sample Latitude for the 3500 km s sample 

0.015 - 



| 0.010 



L = 248° 




0.015 - 



B = 5' 



;o +36° 



! 0.010 




Degrees 

Longitude for the 4500 km s" 1 sample Latitude for the 4500 km s" 1 sample 



L = 223° tl 



B = -10' 





0.015- B = -3' 



Fig. 12. — The one dimensional distribution of the polar quadrupole vectors 
for the different samples, with the 68% limits indicated by a dashed line and 
the median value by a dotted line. 
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Fig. 11. — The one dimensional distributions of the dipole vectors for 
the different samples, with the 68% limits indicated by a dashed line and the 
median value by a dotted line. 



Sample v r [km s ] 




/ 




b 


3500 


280 -88 


289° 


+23° 
-25° 


-11° 


+24° 
-17° 


4500 


27Q +57 
Z -79 


285° 


+15° 
-20° 


-10° 


+15° 
-14° 


HF 


? OQ +70 


281° 


+21° 
-24° 


14° 


+16° 
-15° 



Table 1 

Amplitude and direction of the dipole vector for the 
three different samples. 



Sample 


v r [km s ! ] 




/ 


b 


3500 


654 +129 
0:54 -144 


248° 


+49° 
-55° 


co +36° 
3 -45° 


4500 


575 ^ 37 


223° 


+45° 
-44° 


-10° +35 ° 
—35° 


HF 


522 + 127 


236° 


+35° 
-41° 


oo +27° 

3 _ 31 o 



Table 2 

Amplitude of the quadrupole, and direction of the 
polar quadrupole vector for the three different 

SAMPLES. 



|Jha et alJ (l2007[) also provide a crude estimate of the dipole 
amplitude and direction from a subset of 69 supernovae in 
their sample which is roughly compatible with our 4500 sam- 
ple . In the coordinate system of the Local Group they find 
a velocity of 541 ± 75 km s -1 towards a direction of (/,&) = 
(258° ± 18°, 51° ± 12°). If we transform our dipole term from 
the HF sample to the same coordinate system, using the Lo- 



cal Group velocity derived in lRauzy & Gurzadyanl (119981) . we 
find a velocity of 5 16 km s" 1 towards (b,l) = (248°, 51°). Both 
the amplitude and direction are compatible with the JRK value 
at la. Our derived amplitude is slightly lower (although not 
significantly so), because in our fit the quadrupole term ac- 
counts for part of the velocity. 

As noted before, we do not discuss the local monopole 
term By subdividing supernovae into l ow and high redshift 
bins (Jha et al. 20071: IZehavi et al.ll 19981) . a significant varia- 
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tion in the local Hubble parameter has been detec ted. This 
does not have any impact on our resu lts (see section 15. II) , and 
since it was discussed thoroughly in lJha et all (|2007[) we refer 
to that paper for further details. 

6.2. Galaxy surveys 

Results from galaxy velocity surveys on scales of order 
4000-6000 km s" 1 generally a gree that the magnitude of the 
dipole is of order 30 km s -1 in the direction / ~ 300°, b ~ 
20° (see for instance IZaroubil (120021) and references therein). 
This result is compatible with the SN la dipole direction and 
magnitude within 2a. 

A reconstruction of the very local velocity field (< 
3000 km s -1 ) was done by lTonry et alJ (l200Qh measuring sur- 
face brightness fluctuations in 300 early type galaxies, pre- 
dominantly in groups and clusters. They used an explicit flow 
model with a Virgo Attractor and a Great Attractor which con- 
tain the main local mass concentrations. Furthermore they 
added dipole and quadrupole terms to account for the gravita- 
tional pull and shear from large scale structure further away. 
They find a very low value of the dipole (~ 150 km s -1 ) and 
the quadrupole polar vector (~ 50 km s -1 ), but it may be re- 
lated to having the dipole in the same direction as the attrac- 
tors, and the attractors accounting for the major part of the 
shear (quadrupole term) in the model. 

The dipole has also been measured using velocity field re- 
construction of the 2mass catalogue. At a distance of 4000- 
6000 km s -1 the dipole direction is found to be roughly / ~ 
250 °, b ~ 35-40°, again compatible with our re sult within 
2a (lErdogdu et al.ll2006allbl : IPike & Hudsonll2005l) . 

6.3. Clusters 

Cluster samples like SMAC probe larger distances, and 
find directions which are generally compat ible with the SN 
la result. For example Hudson et alJ (120041) find / = 260° ± 
13°, b = 0° ± 11°. However, they find an amplitude of 
687 ± 203 km s -1 , significantly higher than our result (al- 
though again compatible at 2a). 

7. DISCUSSION 

We have analysed mock supernova surveys in order to study 
the number of supernovae needed to probe the large scale ve- 
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locity field of the local universe, quantified in terms of the 
angular power spectra as a function of redshift. We then 
proceeded to use the best available database of low-redshift 
supernovae, the JRK sample, to probe the local dipole and 
quadrupole of the velocity field at three different distances. 
The present method has several advantages over galaxy sur- 
veys. The uncertainty on each individual supernova luminos- 
ity is much smaller than the systematic uncertainties in deter- 
mining galaxy luminosities so that a much smaller sample is 
sufficient. We find that 

• With two different models for the type la supernova 
rate, the resulting mock surveys only differ at the per- 
cent level. Hence, using SN la to probe the underly- 
ing velocity field is robust with respect to assumptions 
about the supernova environment. 

• For the dipole we find a result which is consistent with 
galaxy surveys at the same Hubble flow depths. 

• The quadrupole is comparable in value to the dipole, 
indicative of a significant shear in the local velocity 
field, in accordance with our mock catalogues. It has, 
to our knowledge, not been measured before at these 
distances. 

• With the present sample size of almost 100 supernovae 
the precision of the dipole measurement is comparable 
to that in galaxy surveys using thousands of galaxies. 

Finally, we note that new surveys like Pan-STARRs, 
SkyMapper and LSST will measure about 10,000 type la su- 
pernovae at z < 0.1 per year, and if proper light curves and 
redshifts can be measured for even a small fraction of these 
events they will provide an extremely powerful tool for study- 
ing the dynamics of the local universe. 
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